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ABSTRACT 

Stellar streams provide unique probes of galactic potentials, with the longer streams 
normally providing the cleaner measurements. In this paper, we show an example of 
a short tidal stream that is particularly sensitive to the shape of the Milky Way's 
dark matter halo: the globular cluster tidal stream NGC 5466. This stream has an 
interesting deviation from a smooth orbit at its western edge. We show that such 
a deviation favours an underlying oblate or triaxial halo (irrespective of plausible 
variations in the Milky Way disc properties and the specific halo parametrisation 
chosen); spherical or prolate halo shapes can be excluded at a high confidence level. 
Therefore, more extensive data sets along the NGC 5466 tidal stream promise strong 
constraints on the Milky Way halo shape. 
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1 INTRODUCTION 

Pure dark matter simulations of galactic halos in our cur- 
rent ACDM paradigm (cold dark matter with a cosmologi- 
cal constant) predict triaxial halo shapes (e.g. Dubinski & 
|Carlberg||1991[|Jing fc Suto||2002[ ). However, if baryonic ef- 
fects are included, the shape becomes more axisymmetric 



and aligned with the gas/stellar disc (e.g. Dubinski 1994 
[Debattista et al.||2008"l [Kazantzidis et al.|[2010[ ). Hence, de- 
termining the shape of the Milky Way halo constrains both 
our cosmology, and our current galaxy formation models. 
The shape of the Galactic potential can also be used to con- 
strain alternative gravity models ( [Read fc Moore|[2005[ ). 
Stellar streams are a powerful tool for probing the Milky 



Way (MW) potential (|Lin fc Lynden-Bell||1977 Ibata et al 



2001 
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, is a text- 


book 


04( jJohnston et al. 


2005 


Law et al. 


2005 


IFellhauer et al.|2006 |Belokurov et al. 


2006 


Law et al. 


2009 


Law fc Majewski 


2010 


1, but its large width, its bi- 



furcation, and the unknown properties of the progenitor 
galaxy introduce large systematic errors to the recovered 



MW halo mass and shape (Penarrubia et al. 2010 20111 



Thinner streams avoid these complications (Lux et al. 2012, 
in prep), but so far have not yielded constraints on the Milky 
Way halo shape owing to one or more of; close proximity or 
low inclination with respect to the MW disc; short length; or 
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poor data quality (c.f. Newberg et al. (20101 for the Orphan 
stream and [Willett et al.| (|2009|); [Koposov et al.| (|2010|) for 



the GDI stream). 

Using only angular positions and radial velocities of thin 



streams, Varghese et al. (20111 find that the observation of 
turning points (apo-/pericentres) is crucial for their method 
to constraint the shape/mass of the probed potential. An ex- 
ample where a potential low apocentre can provide a strong 
constraint is the globular cluster stream NGC 5466. Inde- 
pendently of each other, two groups reported evidence for 



this rather tenuous tidal stream ( Belokurov et al. 2006 Grill 
mair fc Johnson|2006l ) that is an order of magnitude fainter 
than the tidal stream associated with the globular cluster 



Pal 5 ( jOdenkirchen et al.|2003l [Grillmair fc Dionatos|2006 1 

Both groups used the Sloan Digital Sky Survey (SDSS) data 
but different extraction techniques to identify the stripped 
cluster stars. Belokurov et al. ( 2006 1 report a ~ 4° stream 



using neural networks to extract the probability distribu- 



tion of the cluster stars, while Grillmair fc Johnson (20061 
found a 45° stream using an optimal contrast matched-filter 
technique. Because of its limited data set (on-sky positions 
only) and the disputed length of the stream the NGC 5466 
stream has not been used to constrain the Milky Way halo 
shape so far. 



Fellhauer et al. (20071 have modelled the stream nu- 
merically in a spherical halo potential including a disc and 
a bulge. They predict long (> 100°) and faint tidal tails 
for NGC 5466, confirming the extent and attenuation of the 
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Figure 1. The smoothed, summed, weighted image of the SDSS field for NGC 5466 after subtraction of a low order fit (bacltground data 
taken from |Grillmair fc Jo hnson (2006), Figure 1). The magenta error bars mark the position of the stream that we consider as clearly 
connected with the globular cluster. The white square indicates the location of NGC 5466, while the white diamond shows the location 
of NGC 5272. We only trace the NGC 5466 stream for a > 192°, after which we consider its location to be ambiguous. 



stream claimed by Grillmair & Johnson (20061. However, 



they are unable to reproduce the puzzling inconsistency with 



a smooth orbit as noted by Grillmair Sz Johnson (2006j for 



the area observed by SDSS at a = 224° and is considered 
ambiguous westward of a = 192°. This ambiguous part of 
the stream corresponds to the deviation from a smooth orbit 



the west end of the stream (a < 190°). Grillmair & Johnson as noted by Grillmair & Johnson (20061 and is the subject 



(20061 speculate that this could owe to either irregularities 



in the Galactic potential, a recent encounter with a massive 
object or a confusion between streams. However, both the 
an alysis by|Fellhauer et aL](j2007) and the orbits calculated 
by Grillmair & Johnson ( 2006 1 assume a spherical shape for 



the Milky Way halo potential. In this work, we show that 
the deviation from a smooth orbit can be explained by a 
non-spherical/non-prolate halo potential. 

This letter is organised as follows: In section !|2]we de- 
scribe the existing data and in !|3]our method. The results 
are presented in ij4]the results and discussed in i|5] 



2 DATA 



For our analysis, we use the data provided by |Grillmair fc| 
Johnson ( 2006 1 as they report longer tails for NGC 5466 
than [Belokurov et al.| ( |2006[ ). Figure [T] shows the data from 
their work with coloured error bars overlaid to indicate the 
positions of the stream as employed in this worlQ These 
positions range from a £ [192°, 224°] as the stream leaves 



^ Note that the true extend of the stream is disputed. However, 
we assume in this work that the stream truly has the extend 
indicated in Figure [l] 



of our investigation. 

For the initial conditions of our orbit integration, we 
use the on-sky position (indicated by the white square in 
Fig. [T]), distance and radial velocity for the globular clus- 
ter as given in the Harris catalogue (Harris||1996 1, i.e. a — 



14 05 27.29 hms, 5 = 28 32 04.0 hms, d = 16.0 ± 1.6 kpc, 
Vr = 110.7km/s. The Galactic centre distances is then 16.3 
kpc. The proper motion data have been taken from ( Grill- 
mair fc Johnson||2006[ ), 



I.e. [u,v,wj 



[290, -240, 225] km/s. 



We assume ±3mas/yr errors for the proper motions (Munn 
et al.|2004 K Both the distances and the proper motions have 
relatively large errors in comparison to the angular position 
and radial velocities. They are varied within 3a and 2a er- 
rors, respectively. Note, that the orbital plane of NGC 5466 
is roughly aligned with the Sun-Galactic centre axis and per- 
pendicular to the disc, i.e. we observe the cluster from nearly 
within its orbital plane. 



3 METHOD 

We analyse possible orbits of the NGC 5466 globular cluster 
that are consistent with the angular position of its tidal tails. 
Specifically, we investigate what effect the halo shape has on 
the angular positions of possible NGC 5466 orbits. 

The orbits of NGC 5466 are derived by integrating test 
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Figure 2. The plots show the orbits of accepted models in the MCMC chain for the real NGC 5466 data and different halo shapes: 
Spherical (qi = = 1), oblate (51 = 1, ^ ^ 0.95), prolate (gi = 1, 1.05 ^ s£ 1.8), triaxial (1.05 ^ 91 ^ 1.8, ^ s? 1.8). 

Orbits marked by black dashed lines have an apocentre in the Galactocentric coordinates, that appears within the field of view in 
projection. Note, that nearly all orbits with some deviation from a smooth orbit have an apocentre in projection within the field of 
view. Orbits without such a low apocentre are marked with grey lines. Note that we have omitted 90% of these for clarity. The relative 
percentage of orbits with and without a low apocentre is not a converged quantity for the oblate and triaxial halo potentials and can 
vary widely depending on the individual MCMC chain (percentages > 50% are possible). However, we found no such 'kink' orbits for 
the prolate halo potential, and just of order one out for 4 MCMC chains with 10^ samples for the spherical halo potential. (Note that 
the spherical halo potential model is still slightly oblate due to the presence of the disc.) 



particles from the position of the globular cluster forward 
and backwards in time. The initial conditions are taken from 
the globular cluster data as described in ^ keeping the an- 
gular positions and radial velocities constant while varying 
the distance and proper motion data within their errors. We 
integrate the test particle orbits using the Drbit_Int code 



described in Lux et al. (20101. 



We transform the Galactocentric coordinates into ob- 
servable coordinates using M. Metz's tool bap.coord^ 
i Metz et al.|2007 l adopting 8.0 kpc as the sun-Galactic cen- 
tre distance and a circular velocity of 220km/s at the solar 
position. 

For our Milky Way potential model, we fo llow |Law fc| 
Majewski| ( |2010[ ) and use a |Miyamoto &: Nagai| < |1975[ l disc 

G A/disc 

R'^ + {a + Vz^ + b'^y' 



Hernquist ( 1990 1 bulge 



GMbu 



bulge 



Igc 



r -I- c 

and a cored triaxial logarithmic potential 



(2) 



"l>haio = -Uhaio ln(Cia:: + C2y + C-jxy + (z/q:,) -|- r'haio) (3) 



where the constants are defined as 



C2 = 



+ 



+ 



1\ 



Ca = 2sin0cos(^ 



(4) 
(5) 
(6) 



The disc mass is set to Maisc = 1-0 x 10 M© with a scale 
length a = 6.5 kpc and scale height h = 0.26 kpc. The mass 



|http : / /vwii . astro ■ unl-bonn . de /~mmetz/py/docs/mkj _libs/| 
public/bap. coords-module .html 



of the bulge is fixed at Afbuigc = 3.4 x 10^ Mq and its scale 
length c — 0.7 kpc. The logarithmic potential is chosen in 
such a way, that the circular velocity of the total potential 
at the position of the sun (Rq = 8 kpc) is equal to 220 km/s. 
The core of the logarithmic potential is set by rhaio = 12 kpc 
and the shape parameter g2 = 1 is kept fixed. By changing 
gi G [1,1.8], (j) G [0,7r] and qz € [1,1.8] we can cover the 
whole range from purely spherical to fully triaxial halo po- 
tentials. Our results are not sensitive to this choice of po- 
tential functions and parameters (see discussion in section 
|4f. 

We analyse orbits in spherical (gi = Qz ~ 1), oblate 
(gi = 1, 072 s; gz s; 0.95), prolate (gi = 1, 1.05 ^ g^ ^ 
1.8) and triaxial shaped halo potentials (1.05 ^ gi ^ 1.8, 
•\/l/2 ^ g^ ^ 1.8). Note, that both the oblate and prolate 
halo potential have their symmetry axis aligned with the 
symmetry axis of the disc. All other configurations are in- 
cluded under the label "triaxial". In each potential we cal- 
culate orbits consistent with the tidal stream positions as 
marked in Figure I^We use a Markov Chain Monte Carlo 
technique (MCMCf^J with a varying step-size to effectively 
scan the parameter space varying the proper motions and 
distances of the initial conditions as well as the shape pa- 
rameters, if appropriate. 

Note that only specific globular clusters streams with 
negligible stream-orbit-offsets can be analysed with the sim- 
ple test particle fitting technique (Lux et al. 2012, in prep). 
Due to the special alignment of the NGC orbital plane with 
our position, we expect the stream offset from the orbit to 
be mainly in the distances and not in the angular positions. 
Then the expected stream-orbit-offset for this globular clus- 
ter stream is significantly smaller than the current distance 
measurement errors. However, here we do not aim to quan- 
titatively constrain the individual shape parameters of the 



^ The MCMC method is a probability distribution space sam- 
pling method employing a random walk to efficiently scan com- 
plex parameter spaces. The exact details of our implementation 
are described in Lux et al. (2012, in prep). 
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Figure 3. The plots show the radial velocity data and helio- 
centric distances (assuming a Sun - Galactic centre distances of 
8kpc) of the accepted models in the MCMC chain for the oblate 
halo potential. The on-sky positions for these orbits are shown 
in Figure |2] Orbits marked by black dashed lines have an apoc- 
entre in the Galactocentric coordinates that appears within the 
field of view in projection. Note that nearly all orbits with some 
deviation from a smooth orbit in the on-sky positions have an 
apocentre in projection within the field of view. Orbits without 
such a low apocentre are marked with grey lines; we have omitted 
90% of these for clarity. The two dashed vertical lines mark the 
range of the NGC data as employed in this work. 



Milky Way halo, but simply point out the different types of 
orbits possible in various halo shapes. This should not be af- 
fected even if offsets in the angular positions occur. We test 
the convergence of our MCMC chains by running at least 
4 chains with 10^ iterations each starting from different pa- 
rameter sets. 



4 RESULTS 

Figure [2] shows the orbits consistent with the NGC 5466 
stream data for the four different halo shapes: spherical 
(171 = — 1), oblate (gi = f, \/l/2 ^ qz ^ 0.95), pro- 
late {qi = 1, f .05 < gz < 1.8) and triaxial (1.05 ^ gi < 1.8, 



•\/l/2 ^ ^ 1.8). We find that, while the spherical and 
prolate halo shapes only allow smooth orbits, in the oblate 
and triaxial halo shapes 'kinky' orbits occur that can eas- 
ily explain any deviations from a smooth orbit as seen by 



Grillmair & Johnson ( 2006 1 at this end of the stream. Unfor- 
tunately, the deviations can either curve 'upward' or 'down- 
ward' making a prediction for the further location of the 
stream without any further data difficult. In both Figures 
[2] and |3] orbits with a low apocentre have been marked by 
black dashed lines. They show that these orbits not only 
distinguish themselves from the 'smooth' orbits in angular 
positions but also in the distances and radial velocities along 
the stream. 

To test the robustness of our results, we repeated our 
analysis using a triaxial NFW potential (parametrisation as 
in Guedes et al. 2009 1 instead of the logarithmic potential 
for our halo (omitted for brevity). As for the logarithmic 
halo potential, only oblate and triaxial halo shapes allow 
the existence of kinky orbits. 

The existence of the disc potential already introduces a 
certain oblateness in the overall potential even for spherical 
halo shape. Therefore, one might suspect a degeneracy be- 
tween disc and halo potential and expect to see 'kinky' orbits 
even for fully spherical halos. However, we use a conserva- 



tively high disc mass ( Binney & Tremaine 2008 1 and find 



the occurrence of a 'kinky' orbit for a spherical halo shape 
negligibl^ For the prolate halo shape, 'kinky' orbits have 
never been found in our analysis. Therefore, the oblateness 
of the disc potential is not sufficient to mask the halo shape. 
Additionally, we tried replacing the Miyamoto-Nagai disc 
with an exponential disc to test whether this more realistic 
disc potential influences the oblateness of the total potential. 
We find that this is not the case and the amount of 'kinky' 
orbits for spherical halos - given our assumed data - is as 
negligible as for the Miyamoto-Nagai disc (~ 1 occurrence 
in 4 MCMC chains with 10^ samples). 

Note, that the parameter space searched by the MCMC 
chains is degenerate due to the limited data set we are try- 
ing to fit. Therefore, the amount and the type of 'kinky' 
orbits found for the oblate and triaxial halo potentials can 
vary significantly between individual chains. However, for 
the prolate halo potentials they are robustly consistent with 
zero, while for the spherical halo case such orbits are ex- 
tremely rare. For the oblate and triaxial plots in Figure [2] 
we decided to show MCMC chains with rather few 'kinky' 
orbits as the less crowded plots are more instructive. Un- 
fortunately, the nature of the 'kink' (even combined with 
distance or velocity information along the stream) does not 
allow a distinction between the oblate and triaxial halo po- 
tentials. 

We find that the low apocentres are 8-15 kpc away from 
the globular cluster (3D distance). Montuori et al. (20071 



find that GCs on eccentric orbits near apocentre are ex- 
pected to have multiple tails around the globular cluster. 
This could create further deviations from a smooth stream. 
Additionally, streams at apocentre are compressed in both 
angular positions and radial velocities and therefore might 
yield poorer constraints on the underlying potential than 



Please note, that this might depend on the disputed extend of 
the stream and the specific data points chosen for this work. 
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stream further away from apocentre ( Eyre fc Binney|2009 1. 
However, a lower apocentre also means a lower eccentricity 



and Johnston et al. ( 2008 1 find a more stream-like morphol- 



ogy for less eccentric orbits. As low apocentres only exist in 
oblate/triaxial halo potentials, this does not effect our re- 
sults. Testing the effects of a close apocentre on the stream 
morphology with N-body models is therefore beyond the 
scope of this paper and will be postponed to future work. 



5 DISCUSSION AND CONCLUSIONS 

We investigate a possible explanation of the apparent devia- 
tion from a smooth orbit of the globular cluster tidal stream 
NGC 5466 westward of a > 192°. We integrate orbits in 
a variety of halo shapes consistent with the known loca- 
tion of the stream and find that only for either oblate halo 
potentials (with respect to the disc) or fully triaxial halo 
potentials such deviations are possible. If this deviation is 
verified, it places a strong constraint on the shape of our 
Milky Way halo and both spherical and prolate halo shapes 
can be excluded with high probability. 

The observed deviations are highly correlated with low 
apocentres in the globular cluster orbits and can curve either 
way from the current orbit. This makes a prediction for the 
halo shape without additional data difficult. However, these 
low apocentre orbits not only show a very distinct behaviour 
in on-sky positions, but also in radial velocities/distances 
along the stream. Therefore, further maps of the westward 
end of the stream as well as radial velocities/distances along 
the whole known stream are the key to unraveling the Milky 
Way halo shape using the NGC 5466 stream. 
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